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Abstract 



■ Cylindrically-symmetric solutions in Conformal Gravity are investigated and several new 

ov 

\Q . eralizations of the Melvin solution and cosmic strings of the Abelian Higgs model. The 



solutions are presented and discussed. Among them, a family of vacuum solutions, gen- 
eralizations of the Melvin solution and cosmic strings of the Abelian Higgs model. The 
Melvin-like solutions have finite energy per unit length, while the string-like solutions do not. 
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1 Introduction 

Conformal Gravity (CG) [1] was proposed as a possible alternative to Einstein gravity ("GR"), 
which may supply the proper framework for a solution to some of the most annoying problems of 
theoretical physics like those of the cosmological constant, the dark matter and the dark energy. 

The gravitational field in CG is still minimally coupled to matter, but the dynamical basis 
is different: it is obtained by replacing the Einstein-Hilbert action with the Weyl action based 
on the Weyl (or conformal) tensor C K \^ U defined as the totally traceless part of the Riemann 
tensor: 
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so the gravitational Lagrangian is 

where a is a dimensionless parameter. The gravitational field equations take the following form: 

WW = f ^W (1-3) 
where Ij^ is the energy-momentum tensor and is the Bach tensor given by: 

WW = ^V M V^ - V A V A i?^ + ^(i? 2 + V A V A i? - m KX R KX )g^ + 2R KX R I1KI/X - ^RR^. (1.4) 

It was suggested (see [1] and references therein) that while CG agrees with Newtonian grav- 
ity in Solar System scales, it further produces a linearly growing potential that could explain 
galactic rotation curves without invoking dark matter. It was further argued that accelerating 
cosmological solutions of CG describe naturally the accelerated expansion of the universe thus 
removing the need for dark energy. 

On the other hand, CG has been criticized from several aspects both phenomenological and 
formal. Arguments in favor of the need of dark matter come from observations of the unusual 
object called "bullet cluster" [2, 3] whose dynamics seems very difficult to understand without 
assuming a weekly interacting dark component. 

More specifically, several authors claim that predictions in the weak field limit of CG disagree 
with solar system observations [4] , yield wrong light deflection [5] and that the exterior solutions 
cannot be matched to any source with a "reasonable" mass distribution [6]. Other authors find 
evidence for tachyons or ghosts [7] or raise the fact that only null geodesies are physically 
meaningful in this theory since the "standard" point particle Lagrangian is not conformally- 
invariant [8]. 

Counter arguments to some of these objections were also published [9, 10], and the matter 
is, to our view, still waiting for a consensus. 

A somewhat different approach to the subject, was suggested recently [11] and solves natu- 
rally some of the above-mentioned problems, like that of the meaning of time-like geodesies. It 
is based on extending CG into a scalar-tensor theory by introducing an additional real scalar 
field, so the conformally-invariant point particle Lagrangian will be 

L pp = -S^Jg^x^x" (1.5) 

where S is a real scalar field with the usual conformal transformation laws. 

It is therefore very much required to investigate further the predictions and consequences 
of CG in its purely tensorial formulation as well as in its scalar tensor extension as much as 
possible. 
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Following our previous studies of spherically-symmetric solutions [11] including the non- 
Abelian case [12] in CG, we move now to cylindrically-symmetric solutions and especially to 
cosmic strings and we start an investigation of their properties in this context. The present 
report may be considered therefore as a third step of this program. 

Cosmic strings in CG will be the subject of the last section (Sect. 5) of this paper. After 
a short general discussion of the field equations of the self-gravitating Abelian Higgs model in 
the cylindrically symmetric case (Sect. 2), we will start (Sect. 3) with vacuum solutions of 
CG and line source solutions. Then we will discuss purely magnetic solutions (Sect. 4), i.e. 
the analogs of the Melvin solution of ordinary GR and finally we will get to cosmic strings of 
the Abelian Higgs model. We will see that unlike their GR counterparts, the Abelian Higgs 
string-like solutions in CG although localized, do not have finite energy per unit length. 

Since a large part of our discussion will be based on the Abelian Higgs model, we end this 
section with a few words about its coupling to CG. 

Among all the higher order gravitational theories, CG is unique in the sense it is based 
on an additional symmetry principle. The conformal symmetry imposes severe limitations on 
the allowed matter Lagrangian, but the Abelian Higgs model is essentially still consistent with 
the conformal symmetry provided the scalar field "mass term" is replaced with the appropriate 
"conformal coupling" term which introduces a coupling to the Ricci scalar R. The matter 
Lagrangian is therefore 

C m = \(D^y(D^) - li?|$| 2 - - \f^\ (1.6) 

and the resulting field equations are 

+ \\<S>\ 2 $ + — $ = (1.7) 

1 P 

The energy-momentum tensor T^ u is given by: 

Tpu = Tft inima » + l - (^V A V A |<I>| 2 - V M V„|<I>| 2 - GV|<D| 2 ) (1.8) 

rpijninimal) k e i n g ordinary ("minimal") energy-momentum tensor and G^ v is the Einstein 
tensor. 

We can see already at this stage that in order to obtain string-like solutions with |<E>| ap- 
proaching asymptotically a constant, the geometry cannot be asymptotically flat. The most that 
can be obtained are solutions which have asymptotically constant and negative Ricci scalar. Ac- 
tually, a constant Ricci scalar is not a "gauge invariant" concept in CG; neither is a constant 
scalar field. Thus, it is only a matter of convenience which can be obtained by a proper gauge 
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choice. Moreover, we would like the gauge choice to be consistent with symmetric vacuum 
solutions. Thus, we will choose a gauge which allows to impose the condition 

-»• ~S£ for r -»• oo (1.9) 

where k is a positive parameter. We will also limit our solutions to those exhibiting boost 
symmetry in one direction (say, z), as for ordinary cosmic strings. These restrictions will simplify 
considerably the very cumbersome expressions of the components of the Bach tensor and will 
enable a clear physical picture. 

The scalar-tensor extension of CG is done along these lines. The gravitational Lagrangian 
(1.2) will be modified to 

C g = \ {-\c K ^ v C^ + ^V X SV X S - ±RS 2 - V -S^j (1.10) 

where S is a real scalar field with the usual conformal transformation laws and v is a possible 
self-coupling parameter. This should go together with postulating the conformally-invariant 
point particle Lagrangian, (1.5). Point particles couples in this picture to the "physical metric" 

However, when the Higgs field is present, we may identify S with |$| so having a more 
economic model where the Higgs field is responsible for the mass of point particles also in this 
classical framework. 



2 Cylindrically-Symmetric Solutions 



The general cylindrically symmetric line-element has the form: 

ds 1 = B 2 (r)dt 2 - M 2 {r)dr 2 - L 2 (r)d^ 2 - K 2 (r)dz 2 

and the general expression for the Ricci tensor turns out to be 
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Since all the solutions that we will consider in this work can be thought technically as 
special cases of the solutions of the Abelian Higgs model coupled to CG, we discuss now the 
field equations of this system. 
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The general cylindrically symmetric matter fields are parametrized by 



$ = /(r)e imv ; cA^dx* = A(r)d<p. 



(2.4) 



Using the above parametrizations, the field equations for the scalar and vector fields get the 
following simple form: 
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(2.5) 
(2.6) 



The gravitational part of the field equations is less simple due to the complicated form of the 
Bach tensor. First we write explicitly the non-vanishing components of Ty and in order to 
simplify matters we define the following partial energy densities: 



iff _ 1 (A!_ 

£s ~ 2 ' £v ~ 2 \L 

In terms of these the components of T„ will be 



1 (m-A 

2 V L 



u 



A 4 
4 J ' 



(2.7) 



1 



rpz 
1 Z 



1, 



— £u £"s H~ £s 



u 



f 2 



RL 



3 V B L 



6 V * 6 



+ 



3 L 



(2.8) 



where the Ricci tensor components are given in Eq. (2.2). As an easy check one can verify that 
indeed the sum of these components vanishes as it should. Note that unlike in GR, this does 
not force a vanishing Ricci scalar since the gravitational field equations are given by Eq. (1.3). 
In the present case they reduce to two components of (1.3) which we may choose to be the tt 
and rr components: 
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(2.9) 



We still refrain from writing explicitly the components of W„ anticipating further simplifications. 

In order to find solutions for this system, we have to fix the arbitrariness of the radial 
coordinate and the arbitrary rescaling of the metric due to the conformal symmetry. This, 
together with the boost symmetry in the tz plane leaves one independent metric component. So 



only one of the gravitational field equations has to be solved. We chose to solve the lower order 
rr equation of (2.9), and used the higher order tt equation as a consistency check. 

The choice of the conformal rescaling further has to be compatible with the condition oc 5^ 
asymptotically (see (1.9)). This is not obvious even in vacuum: for instance, in the cylindrical 
version of the "Mannheim gauge" [1] M(r) = 1/B(r) ,L(r) = r with K(r) = B(r), the vacuum 
equations are inconsistent; in the gauge M(r) = l,K(r) = 1 , the vacuum equations are consis- 
tent but the solutions are such that the diagonal components of Ricci tensor depend linearly on 
r for r — > oo. 

One convenient ansatz which solves these problems is 

M = 1 , L = — , K(r) = B(r) (2.10) 

where we also resort to dimensionless coordinates. 

Moreover, all the thin string (line source) solutions of GR with a negative cosmological 
constant [13] which solve for r > 

K = -\% (2-11) 

satisfy also the CG vacuum equations = 0. Among those solutions is also found the (four 
dimensional) AdS soliton [14] (see also [15, 16, 17]) which is a cylindrically-symmetric regular 
solution of the same equation (2.11) and is therefore distinct from AdS (anti-de Sitter) space. 
In this gauge the components of the Bach tensor are 
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while the fourth one, W$ can be obtained immediately from the identity Wjt = 0. 
In terms of B(r), the Ricci scalar and tensor take the form 
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Returning now to the general Abelian Higgs system we are therefore left with the rr com- 
ponent of equations (2.9): 
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2 V ^ / 2 2 \ L y 4 6 V r 2 / 3 V 5 L 

which should be solved together with the two equations (2.5)-(2.6). 



3 Vacuum Solutions 

We start by considering the vacuum equation W^ v = associated with the metric (2.1) in the 
gauge (2.10). The relevant equation for the function B(r) is just a special case of Eq (2.15), 
leading to 



2gg) 2B^B" 2B"'B" 2 2B'" B" 2B'"B' B'" 2 AB" 2 AB' 2 B" B' 4 _ 

~3B + 3B^ BB^ + 3B 2 + W* 2 + ~W 2 ~ 35 3 + W 4 ~ ( ^ 

Equation (3.1) has very interesting properties, namely it is autonomous and invariant under 
rescaling of B and of r. 

The above-mentioned AdS soliton is just one member of a continuous family of solutions of 
(3.1) (as well as of (2.11)) given by 

B(r) = (acosh(kr) + 6sinh(/cr)) 2/3 . (3.2) 

The case k = \/3/2 and 6 = corresponds to the AdS soliton which satisfies Eq. (2.11) with 
k = 8. The other values of k with 6 = correspond to thin string solutions which are the 
counterparts of the conic solutions for k = 0. The solutions with a = are analogous to the 
"Melvin branch" of thin strings. These have a power law behavior near the string axis with the 
same powers of (2/3, —1/3, 2/3) as for k = 0. Actually these exact power law k = solutions 
also solve Eq. (3.1) and may be viewed as a limit of (3.2) where k — > and kb remains finite: 

B(r) = (a + 6r) 2/3 . (3.3) 

The complexified version of (3.2), namely 

B(r) = (acos(kr) + 6sin(fcr)) 2/3 (3.4) 



are closed solutions in analogy with the "inverted cones" of the k = case. 
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Since the gravitational field equations of CG are of higher order, it is expected that more 
vacuum solutions exist, i.e. new vacuum solutions which are special to CG. Indeed the following 
family also solves Eq. (3.1) - but satisfies (2.11) only asymptotically: 

B{r) = (acosh(kr) + &sinh(£x)) 2 . (3.5) 

Eq. (3.1) is solved also by their complexified counterparts, 

B(r) = (acos(kr) + 6sin(fcr)) 2 (3-6) 

as well as by the k — > limit solutions whose Ricci scalar vanishes asymptotically: 

B(r) = (a + br f . (3.7) 

We have again a similar pattern of two families of string-like solutions: an open one and a 
closed one, each of them with the free parameters a, b and k. The family of open solutions with 
6 = contains as before a soliton-like regular solution for the specific value k = l/\/2. This one 
is however a new solution special to CG. We write both solitons explicitly: 

K = 8 : B(r) = cosh 2/3 (y/3/2 r) , R(r) = -8 

k = 24 : B{r) = cosh 2 (r/Vtj , R(r) = -12 (l + tanh 2 (r/v 7 ^)) . (3.8) 

It is likely that equation (3.1) is integrable by quadrature but, so far, we were unable to find 
the general solution analytically. One additional explicit solution that we found is B(r) = e ±fcr 
which is actually conformally flat and thus uninteresting for us. In absence of knowledge about 
the general set of analytical solutions we investigated the behavior of the solutions of Eq. (3.1) 
numerically. Exploiting, the different scale invariances of (3.1), the boundary conditions of 
regular solutions can be specified according to 

B(0) = 1 , B'(0) = , B"(0) = 1 , R(t — y oo) = — k (3.9) 

where k is a positive constant. The two regular solutions (3.8) already obey the boundary 
conditions (3.9). Note the residual scaling symmetry of any solution obeying (3.9) realized by 
the one-parameter family of solutions 

B (r) — > B(r) = AB(-^=) , k^k=j. (3.10) 

The combination <5(r) = {BB" — (B 1 ) 2 ) / B has some relevance in the analysis of the solutions. 
In particular the zeros of this function constitute singular points of Eq. (3.1). The occurrence 
of such points renders the solutions problematic. For k < 8, our numerical analysis reveals that 




Figure 1: The two analytical regular solutions of the vacuum equations. Profiles of B(r), L(r) = B'(r) 
and of R(r), S(r). The line R(r) = —8 for n = 8 is not shown. 

the function 5(r) develops nodes and suggests that no regular solution on r G [0, oo] can be 
constructed. The boundary conditions lead to 5(0) = 1. 

The functions B(r) of (3.8) and their derivatives B'{r) = L(r) are plotted on Fig. 1. We 
also added plots of R(r) and 8(r). 

We then attempted to construct numerically the solutions obeying (3.1)+(3.9) for different 
values of k. It turns out that regular solutions exist only for k > 8. The interpretation of 
the numerical solutions relies on the understanding of their asymptotic behavior. The solutions 
obtained obey the following asymptotic behavior 

log(B(r ->• oo)) = 70 + 7ir + 72e- /3r + O(e- 2 ^) , 7 i = ^| (3.11) 

where 70, 72 and /3 are positive constants. The relevant component of the Bach tensor is 

^;(r^oo) = (3 7 i-/3)(7i-/3)^e- 2 ^ , (3.12) 

and, for the Ricci scalar, we have 

R(r -> 00) = -12 7l 2 + 2^(371 - 0X471 - /3)e" /3r (3.13) 

7i 

Actually, this asymptotic form is general enough to be useful for understanding the solutions 
where matter is present (see next sections). 



The numerical results give j3 = 71 for the vacuum solutions which is indeed consistent with 
the vacuum condition Wuv = 0. The quantitative values of the other parameters can be extracted 
from the numerical solutions. The value 72 is related to the parameter 5 = 6(00): 72 = 12^e _7 °. 

Several parameters characterizing the vacuum solutions are presented as functions of k on 
Fig. 2; namely the value B""(0), the difference A = R(0) — R(oo) and the parameter 5 defined 
above. In fact, regular solutions exist for k > 8 but the picture is limited to k < 24. We notice 
that the quantities presented on Fig. 2 all tend to zero in the limit k — > 8 corresponding to the 
AdS soliton. 

Finally, we consider the vacuum solutions of the scalar-tensor CG, namely those which are 
derived from the Lagrangian (1.10). Technically, they are obtained easily from Eq. (2.15) by 
substituting A(r) = and m = and replacing f(r) by S{r)/yfa and A by av. 

We find regular soliton-like solutions by imposing the boundary conditions (3.9) together 
with 

5'(0) = , S{r -»• 00) = S (3.14) 

where Sq fixes the asymptotic value of the Ricci scalar by R(oo) = — 6^Sq. Fig. 3 contains the 
profiles of S(r) and the metric components of a typical solution as well as the "physical metric" 
S 2 9fj,i/ which is the one which directly couples to point particles - see (1.5). We notice that since 
the scalar field S(r) does not change very much, the "physical metric" is quite similar to the 
"conformal metric". Moreover, even the difference with respect to the vacuum metric of the 




K 



Figure 2: The value B""(0), S = 6(00) and A = 1 - R(0)/R(oo) as function of k. Note: k > 8. 
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Figure 3: A vacuum solution of scalar-tensor CG. Profiles of the scalar field S(r), the metric components 
B(r) and L(r) = B'(r) and those of the "physical metric" S(r)B(r) and S(r)B'(r). k = 24 and v = 1 . 



purely tensor theory shown in Fig. 1 is not very pronounced. 



4 Melvin-like Solutions 

Setting $ = in (1.6), we obtain the Weyl-Maxwell Lagrangian and the corresponding cylin- 
drically symmetric equations are obtained by setting f(r) = in Eqs. (2.6) and (2.15). The 
generalization of the Melvin solution to CG can be looked for. In fact the Maxwell equation can 
be integrated directly, leading to A(r) = Q m /B{r) where Q m is an integration constant which 
encodes the magnitude of the magnetic field. Actually, the z-component of the magnetic field, 
given by B z (r) = A' (r) / L(r) , leads to £> z (0) = Q m . Also the choice A — > as r — > oo was taken. 

The different components of the energy momentum tensor can then be computed in terms 
of B(r), namely 

i - l r -e v - 2B4(r) . (4.1) 
The inertial mass per unit length of the conformal Melvin solution can then be computed: 

/f'OO Tjt 
d 2 x^T° = 2irB z (0) 2 J q dr— = KB z {Qf (4-2) 
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4.1 Background Solutions: a = 

Setting first a = 0, leads to solutions in the background of a CG vacuum. The profiles of the 
magnetic functions A, A' corresponding to B z (0) = 1 are shown in Fig. 4 for the cases k = 8 and 
k = 24 corresponding to the two analytic solutions (3.8). The energy density Tq is supplemented 
on the figure. Similar solutions exist for generic values of k, the profiles of the magnetic function 
A depend only weakly on k. 

4.2 Gravitating Solutions: a > 

Next, we investigate the back-reaction of the electromagnetic field on gravity by solving the 
coupled equations (2.6)+(2.15), for a > 0. Using the Maxwell equation to eliminate A, we 
obtain a single equation for B{r) which reads 

2^(4) 2B^B" 2B"'B" 2 2B"'B" 2B'"B' B'" 2 4B" 2 4B' 2 B" B' A aB z (0) 2 

~~3B W' 2 + W' 3 BB' + W 2 + 3B l2+ 3B 2 W 3 + 3B 1 ~ 45 4 ' 

(4.3) 

The constant B z {0) can be clearly absorbed here in a so from now on, we set B z (0) = 1 and 
study the solutions of (4.3) for different values of a and k. 

It should also be noticed that the function B{r) keeps asymptotically the form (3.11) because 
W£ in the left hand side in Eq. (4.3) is asymptotically of order exp(— 2/3r) with /3 = y n/12 




12 3 4 



Figure 4: Profiles of the magnetic potential A, A' and of the energy density Tq in the background of 
the two analytic solutions (3.8). 
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Figure 5: Melvin-like solutions for k = 24: Form of R(r) and S(r) for four values of a. For a < see 
Sect. 4.3. 



while T T r in the right hand side is asymptotically of order exp(— 4/3r). The deviation due to the 
matter field therefore comes out as a next to leading order correction. 

In fact, the deviation of the metric and matter functions with respect to the case a = is 
not substantial. The deformation of the geometry by the matter fields appears more clearly on 
the Ricci scalar R{r) and on the function 8{r) defined in Sect. 3. Fig. 5 shows these functions 
for k = 24 and for three positive values of a. Our results further show that the solutions exist 
up to a maximal value of a, a c ~ 12.1. As another example, for k = 10, we find a c ~ 2.6. The 




Figure 6: Melvin-like solutions for k = 10: Form of R(r) and 6(r) for three values of a. 
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r-dependence of R(r) and 5(r) for k = 10 and for three values of a is illustrated on Fig. 6. This 
feature seems to occur for all values of k and we conclude that for fixed k, the solutions exist 
up to a maximal value of a = a c {n) forming a branch of solutions which we denote 'branch I' 
for later convenience. 

Remarkably, the domain of existence of solutions in the a, k plane is determined by 

a < Oc(k) = 0.675k - 4.1 (4.4) 

with a very good accuracy. It looks as if, for a fixed k, CG cannot support a magnetic field of 
arbitrary large strength. In other words, the value |fiz(0)| has to be smaller than a critical value 
(recall absorbing B z (0) in a). 

In order to complete the pattern of the solutions of Eq. (4.3), we looked for another family 
of solutions. The numerical analysis reveals, indeed, that a second branch (to which we refer 
as 'branch II') of self-gravitating solutions exist for a g]0, a c ]. In the limit a — > a c the two 
branches join, forming a cusp and no solutions exist for a > a c . 

A natural point which raises at this stage is understanding the behavior of the solutions 
of branch II for a — > 0. It turns out that for sufficiently large values of a, the value B""(0) 
becomes negative and the function 5(r) develops a local minimum, say for 8(ro) = S m . The 
results further reveal that, in the limit a — > the solutions of branch II are such that 5 m — > 0. 
As a consequence, the equation develops a singular point in this limit and regular solutions no 
longer exist. 
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Figure 8: Comparison of R(r), S(r) and of the magnetic field B m (r) for the two solutions with n = 12, 
a = 2. The red (resp. black) curves refer to branch I (resp. branch II). 

These phenomena are illustrated in Figs. 7a and 7b where several parameters characterizing 
the solutions are presented as functions of a for k = 12. 

A comparison between the two solutions corresponding to n = 12, a = 2 is shown on Fig. 8. 
The figure reveals that the rather tiny difference occurring for the matter fields has a significant 
influence on the geometry; this appears clearly on the shape of the Ricci scalar R(r). 

4.3 Gravitating Solutions: a < 

The solutions of Sect. 4.1 can be deformed for negative values of a, leading to a continuation 
of the 'branch I' of Sect. 4.2 for a < 0. Indeed, a negative value of a is a "wrong sign" choice 
since it yields a repulsive linear potential of localized solutions in the spherically-symmetric case 
[11]. However, the attractive contribution from the negative cosmological constant is dominant. 
Therefore we do not exclude this possibility. 

The corresponding data is shown on Fig. 7. For fixed k, we observe that the parameter 8 
decreases monotonically with a and reaches 5 = at a minimum value a = a m {K). For example, 
we find a m (12) ~ —1.25 and a m (20) ~ —12.0). For a < a m , Eq. (4.3) develops singular point 
corresponding to 5(r) = and cannot be integrated numerically. The difference of the geometric 
functions 5(r) and R(r) for the two signs of a can be appreciated from Fig. 5. 
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Figure 9: Profiles of a Melvin-like solution in scalar-tensor CG. a = 1, v = 1 and k = 24. 



4.4 Solutions of the Scalar-Tensor Theory 

In order to get the Melvin-like solutions in the scalar-tensor extension of the theory, one may 
start with the general field equations of the Higgs model of sec. 2 - Eqs. (2.5), (2.6) and (2.15). 
We replace A by av and y/af(r) by the real (uncharged) scalar field S(r). Consequently, we 
cross out all the interaction terms containing contributions of the form (m — A)f/L. We solve the 
system numerically and find that the solutions are quite similar to those of the purely tensorial 
one. Fig. 9 presents the gauge potential, the scalar field and the metric components as well as 
those of the "physical metric", S(r)B(r) and S(r)B'(r). Comparison with Figs 1 and 4 shows 
that the additional scalar field does not change much the functions A(r) and B{r). 



5 Conformal Strings 

We finally considered the equations of the Abelian Higgs model coupled to Weyl gravity. Cosmic 
strings [18] are a typical outcome in any field theory which describes matter in the very early 
universe, thus serving very well the purpose of testing the implications of CG. 

Relativistic magnetic flux tubes in flat space were first obtained by Nielsen and Olesen [19] 
and their coupling to Einstein gravity was shown to have important cosmological consequences 
[18]. A complete classification of the solutions of the self-gravitating Abelian Higgs model was 
performed [20] and two branches of solutions were shown to exist: the "string branch" consists 
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of solutions which tend asymptotically to Minkowski space-time with an angular deficit, and 
a "Melvin branch" whose geometry is asymptotically similar to that of the Melvin solution. 
Gravitating cosmic strings in the presence of a cosmological constant, were studied in [15]. 

Cosmic strings in the Abelian Higgs model with a conformal coupling to gravity have been 
also considered before (see e.g. [21] and references therein), but only within the context of GR. 

Now we turn to the analogous solutions within CG, namely solving the full system of the 
three equations (2.5), (2.6) and (2.15). 

Since we are interested in localized solutions, the components of T„ should vanish asymp- 
totically. This, together with the boundary conditions f(r) — > v and A(r) — > m as r — > oo fix 
the asymptotic Ricci tensor to be given by (1.9) with n = 6Xv 2 . Note that unlike the ordinary 
cosmic strings, here the asymptotic value /(oo) = v does not originate from the Lagrangian, 
but is a free parameter which characterizes the various solutions of the given system. 

The system of three equations has to be solved with the boundary conditions 

/(0) = 0, /(oo) = «, A(0) = 0, A(oo) = m (5.1) 

for the matter functions and (3.9) for the metric function B(r). The equations lead to a several 
possible asymptotic forms for the solutions, allowing a few kinds of exponential corrections. 

For \a\ < 2, the form below appears to be the one which best fits the numerical solutions 
obeying the boundary conditions (5.1). For the metric function B(r), we have the same form as 
before : 

log(B(r ->• oo)) = 70 + 7ir + 7 2 e _/3r , 7i = , (5-2) 

where 70, 72 and f3 are constants. The parameter (3 encoding the decay rate of the first correction 
now depends on a. For \a\ < 2 we find a linear dependance 

P = 7i +P « > (5-3) 
where the numerics shows p ~ 0.18 . For the matter fields, the asymptotic form is : 

f(r -» 00) = v(l + F e- 4 ^ r - 1 f i fl 1 l ~ ® e-Pr) , A(r -> 00) = A e~^ sr (5.4) 
V 37i(/3 + 7i) j 

with 71 = \J Xv 2 /2 , s = (1 + \J\ + 8/A)/2 and where Fq,Aq are constants. We limit our 
analysis to the case A = 1, so that s = 2. The third term appearing in / in (5.4) is specific to 
CG because it is related to a source term appearing due to the coupling between the scalar field 
and the Ricci scalar R; it dominates the standard homogeneous piece (the term proportional to 
F ) since \/3\ < 4^. 

The asymptotic form (5.4) leads to several qualitative differences with respect to the Melvin 
case treated in the previous section. The most significant difference resides in the fact that the 
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energy momentum tensor decays according to T® oc exp (— j3r). This implies in particular (see 
(5.3)) 

V=gT$ oc e {2 ^ pa) r . (5.5) 

Within the set of parameters which we were able to explore, the exponent is positive and as a 
consequence, the inertial mass (per unit length) of these cosmic string solutions in CG is diver- 
gent. The presence of the conformal coupling and constant asymptotic curvature ( "cosmological 
constant") has therefore important consequences on the physical characteristics of the cosmic 
strings; similar phenomena were observed in the context of spherically-symmetric topological 
defects, namely for the magnetic monopole [22] and the corresponding CG case [12]. 

Another different feature of cosmic strings with respect to the pure magnetic (Melvin) solu- 
tions appears in the asymptotic behavior of the function 5(r) : 

S(r -> oo) = /3 2 72 e (7 °- p a)r . (5.6) 

This combination of the metric function is not any longer constant for r — > oo. This is observed 
and confirmed by the numerical solutions. 

We first integrated the equations of conformal string numerically in the case a = and 
obtained a family of cosmic string solutions in the background of the vacuum space-time charac- 
terized by the parameter k. The profiles of the different matter functions of the string embedded 




Figure 10: Flux tube solution: Profiles of the matter functions A, A', /, /' and of the energy density 
Tq in the case a = 0, A = 1 and k = 20. 
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in the background of the vacuum solution with k = 20 are presented in Fig. 10 (the energy 
density is also supplemented). The profiles of the matter fields are qualitatively very similar to 
those of the flat case [19] or in the gravitating case [20]. Considering several values of k leads 
to similar plots. 

Finally, we studied the self-gravitating cosmic strings by integrating the field equations for 
a / and found a pattern qualitatively similar to the one obtained for pure electromagnetic 
(Melvin) solutions. The matter fields and metric components of a typical solution are depicted 
in Fig. 11. In this case we let the Higgs field play the role of the scalar S, so we stay with the 
more economic model where the Higgs field itself is related to mass generation already at the 
classical level for point particles. 

Although the profiles of the matter functions deviate a little from the a = case, the 
deformation of the metric function B{r) and some quantities characterizing the geometry are 
more significant. The quantities R(r),5(r) are compared in Fig. 12a for a = and a = ±1. 
The derivatives B", B'" also develop some structure close to the symmetry axis r = 0; this is 
illustrated in Fig. 12b. 

In this case also, the gravitating cosmic string solutions exist only up to a maximal value of 
a. E.g. we find solutions for a < 1.785 (resp. a < 1.2) for k = 20 (resp. k = 16). We strongly 
suspect that, like in the Melvin case, a second branch of solutions exist but it was not attempted 
to construct it. The determination of the domain of existence of the solutions in the k — a plane 
was also left for further investigation. 
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Figure 12: String-like solution for A = 1 and k = 20 and a — (black), a — 1 (red) and a = — 1 (blue): 
(a) Profiles of R(r) and 5(r). Note that the S(r) curve for a — can be fit between the two curves of 
Fig. 1. ; (b) Details of derivatives of B(r). 

6 Conclusions 

In this paper we have investigated cylindrically-symmetric solutions in Conformal Gravity. Ex- 
amining the solutions both by analytical and numerical methods we were able to construct 
several types of solutions which to our knowledge were unknown so far: soliton-like regular vac- 
uum solutions, spaces around line sources, the CG analogs of the Melvin solution and conformal 
self-gravitating flux tubes, i.e. the CG analogs of the cosmic strings 

The analysis was based on the Abelian Higgs model coupled to CG which was studied first in 
several special cases and finally for cosmic strings. The symmetry breaking mechanism which is 
required for flux tube solutions has in the present case a gravitational origin, unlike the ordinary 
case, since conformal symmetry does not allow the usual negative mass term into the Lagrangian. 

The Higgs model was also used as a starting point of extending CG to be a scalar-tensor 
theory which allows a consistent coupling of point particles. 

CG is a rich and interesting theory which may supply answers to some of the most annoying 
problems like those of the cosmological constant, the dark matter and the dark energy. It is 
therefore natural and required to investigate the theory in other domains and search for further 
implications of the theory like in the area of topological defects which are assumed to play an 
important role on structure formation at the early universe. 

As a first step we studied the vacuum solutions of CG which are interesting on their own 
right and also as they give a first indication about the asymptotic behavior of localized self- 
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gravitating structures. We found that the cylindrically-symmetric vacuum open solutions of CG 
have an asymptotical behavior which is similar to that of the AdS soliton. Actually there are 
two kinds of solutions: one with an everywhere constant negative ("AdS-like") Ricci curvature 
(the same as the AdS soliton) and the second, which is a new kind special to CG, has a Ricci 
scalar which approaches a constant only asymptotically. There exist also two kinds of flat or 
asymptotically flat vacuum solutions which we have not studied extensively, since they are not 
compatible with the symmetry breaking mechanism which produces in this context the flux tube 
solutions. 

The second step was to study magnetic solutions in the purely tensorial CG as well as in 
its scalar-tensor extension. We have found Melvin-like solutions that have a finite energy per 
unit length. The solutions separate into two branches and in each of them CG cannot support 
magnetic fluxes with too intense magnetic fields in the core. 

The last step was centered about self-gravitating flux tube solutions in CG. Here the system 
contains from the start a scalar field which may be naturally exploited for consistent coupling 
of point particles, so we did not extend the model by an additional scalar field. We were able 
to solve numerically the field equations and obtained 2 important results: (i) As in the purely 
magnetic solutions, there is an upper bound on the magnetic field in the core above which no 
solutions exist, (ii) Unlike the purely magnetic case, the Ricci scalar ( "cosmological constant") 
influences drastically the asymptotic decay of the scalar fields and makes it impossible to have 
a finite inertial mass. 
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